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C ~ 3 i With the known two loop renormalisation of the magnetic coupling, the 4D results 

of the spatial Wilson loop are compared to the prediction from the magnetostatic 
sector. 

o 

1. Effective theories for hot QCD 
Q-i, Perturbation theory in hot QCD suffer from infrared (IR) divergences. 



These divergences can be best confronted by constructing effective field 
theories for the low-energy dynamics. For the partition function and other 
static observables, the relevant theory is the Dimensionally Reduced (DR) 
■ effective theory 1,2 . 

The main idea of DR is to separate systematically the physical scales 
of the quark-gluon plasma and to construct, for each of these scales, an 
effective action. The procedure is as follows: 

• At distances of order 1/T, physics is described by the non-static 
modes. But, if one wants to go to larger distances, one meets with 
IR divergences, due to massless static modes. 

• At distances of order 1/gT, however, electrostatic modes get 
screened by the Debye mechanism, like in a classical electromag- 
netic plasma. This thermal mass is going to regulate the electric IR 
divergences. At these scale, QCD is described by a 3D electrostatic 
Lagrangian L E that reads: 

L E = Tr(D(A)A Q ) 2 + m 2 E TrA 2 + X E {Tr{A 2 )) 2 + 

+ \ E {Tr{A f \{TrA 2 ) 2 ) + X -TrF 2 + SL E . (I) 



1 
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with as parameters gE(g,T), m,E{g,T), ... and SLe represents 
higher orders operators of relative order <? 4 . 
• Furthermore, at distances of order l/g 2 T, magnetostatic modes 
are the dominant ones. At these distances QCD is described by a 
confining 3D Yang-Mills theory Lm that is non-perturbative. 

L M = ±TrF 2 j +6L M . (2) 

with a gauge coupling 3m(3_e, »)£,...). SLm is of relative order g 3 . 

If one wants to compute the thermal average of an observable, then each 
of these three scales will have its own contribution to the final result. 

From the relative order of the truncation in (1) and (2) we see that gu 
is needed to relative order g 2 . In the next section, I will briefly present the 
result for the coupling gu of the 3D magnetostatic Yang-Mills theory Lm 
at two loop order. 



2. Two-loop determination of the magnetostatic action 

The basic idea behind the effective actions eqns (1) and (2) is that one 
can compute with both in the region of momenta p ~ g 2 T. To know what 
the parameters of the latter are in terms of those of the former requires 
computing two-point functions, three point functions etc. in both theories 
and match them. In the matching the diagrams of the pure 3d Yang-Mills 
theory drop out. 

Here we will follow a well-known shortcut 3 by introducing a background 
field Bi in Le- 

A = —B + Q 

9E 

A = gEQo- (3) 
We calculate the fluctuations around the background in a path integral: 

cxp 1 -S M (B) = f DQ a DQi cxp ( — Se — jTr(DiQi) 2 ) . (4) 

9m J s 

We added a general background gauge term. The resulting action 
Sm(B) is gauge invariant to all loop orders and the renormalization of 
the coupling is identified from the background field two point function at a 
momentum p = 0(g 2 T). 
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cxp (_ J_S M (B)) = cxp ( - J T SM(B))(1 + (F^ + F* r + ...)S M (B)). (5) 
9m 9e 

Fi is the sum of all Feynman diagrams for the two point function of the 
background field with i loops and reads: 



Fi = F* r (Si mP 2 -pip m ). (6) 
This leaves us with the relation, using eq.(5): 

4- = i-*f-*f, (7) 
9m 9e 

with 



ff |jf = __L^ ( 8 ) 
9%F? = - (9) 
This is the main result 4 . 



3. Spatial Wilson Loop 

We want now to test the applicability of 3D physics at medium high T. Let 
us consider an observable which has, unlike the pressure whose dominant 
contribution is the Stefan-Boltzmann term due to hard modes, its dominant 
contribution from the 3D Yang-Mills sector. Such an observable can be the 
spatial Wilson loop in the fundamental representation: 

W{L) = TrPcxp(i <j> gA-df). (10) 

As L is purely spatial, it measures the magnetic flux in the plasma. The case 
where A is in an irreducible representation made of k quarks is described 
in this volume by CP. Korthals Altcs. The thermal average of this spatial 
Wilson loop shows area behaviour with a surface tension cr(T). As it is 
a purely magnetic quantity, we expect from dimensional arguments that 
\F® = c 9m where c is a nonperturbative proportionality constant. Indeed, 
as the average of the loop is due to long distance correlation, hard modes 
will not have any effects on the thermal average. In the same way, for soft 
modes, we can integrate out the Aq field, which is what we have done in 
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the previous section while constructing Lm, as the loop does not depend 
on A ! Finally, we have: 



(W{L)) = exp-crA(L) = 



/ DAW(L) exp-S M (A) 



J DAexp~S M (A) 
and it gives, as SLm is of relative order g 3 , 

a(T) = c 2 g 4 M ^ + 0(g 3 )). 



(11) 



(12) 



Now the aim is to fit with our formula for gu the proportionality con- 
stant c, and to study its extension for finite T. For that, we have to go to 
the lattice. We are going to fit has a function of Indeed, gu is a 
function of gE and m# which are functions of g and T. So, for N = 3: 



19.g 2 



16tt 512tt 2 



with 



^ (Log(^) + Log(^-) + 1.90835...). 



(13) 



(14) 



So thanks to these formulas, one can fit -?= versus i£- in order to deter- 
min c and -r^—. 

We have taken data points for the SU (3) spatial string tension from 6,9 
for T > 2T C . So we have 10 points. The fit is shown in fig 1. 




Figure 1. Plot of -S= versus and fit with our two-loop formula. Each of the 10 
points is shown with its error bar. 



We have found the following results: 
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T^Z = L78 ( 12 ) ( 15 ) 

c = 0.505(11). (16) 

Then we have also fitted data with T > 2.5T C . This leaves us with 6 
points. We obtained 

= 1.57(20) (17) 

1V MS 

c = 0.488(18) (18) 

Let us draw our conclusions. 

Despite the fact that the region of confidence is larger, the value of c is 
still incompatible with refs. 8 ' 7 ' 5 . To see compatibility one obviously needs 
string data for higher T. It might also be recommcndable to use other 
Monte Carlo updates, like the one described in rcf. 11 , which is a version of 
the Luscher-Weisz algorithm 10 . 
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